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Abstract 

The effect of noise on a quantum system can be described by a set of operators obtained from the 
interaction Hamihonian. Recently it has been shown that generalized quantum error correcting 
codes can be derived by studying the algebra of this set of operators. This led to the discovery of 
noiseless subsystems. They are described by a set of operators obtained from the commutant of 
the noise generators. In this paper we derive a general method to compute the structure of this 
commutant in the case of unital noise. 
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I. INTRODUCTION 



Quantum mechanics promises to manipulate information for communication, cryptogra- 



phy and computation in a way fundamentally different from its classical counterpart [21 1 . 
Although it is possible to manipulate small quantum systems in the laboratory, the task to 
do so for large ones is daunting, especially because in absence of control of noise and im- 
perfection of realistic devices the quantum properties of the state are destroyed. Quantum 
error correction methods have recently been discovered which protect quantum information 
against corruption. In particular it was shown that if the error rate is below a certain 
threshold 

mm 

then it is possible implement quantum error correction and obtain 
a fixed probability of success with only a polynomial amount of resources. The work on 
the accuracy threshold assumes certain error models such as independent or quasi inde- 
pendent errors. The theorem shows that under reasonable assumptions, imperfections and 
imprecisions of realistic devices do not create fundamental objections to scalable quantum 
information processing. 

The estimated value of the error threshold is small and, today, an extraordinary challenge 
for experimentalists. From the experimental point of view there is a need to better under- 
stand the details of error models of quantum devices, the form of the Kraus noise operators 
which occur during the evolution of a quantum computation or during the transmission in 
a quantum channel, the correlation between errors, their strength, etc. From the theoretical 
side we need to find efficient ways to estimate the error model, determine the threshold for 
various physically relevant error models and find ways to achieve better threshold. This can 
be done by optimizing quantum error correction procedures or simply finding better ways 
to protect quantum information. 

Recently a unified view of error correction was given which brought together quantum 
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error correcting codes and noiseless (decoherence-free) subspaces 
investigation led to the discovery of noiseless subsystems. A quantum bit, called a qubit, 
is described by a unit vector in a 2-dimensional Hilbert space. Often this 2-dimensional 
space is directly associated with a two-level physical system such as the ground and excited 
states of an atom, a nuclear spin pointing up or down, etc. Noiseless subsystems are instead 
described by a set of operators such as the Pauli matrices, a more general concept than 
a subspace of a Hilbert space Q]. An example where the difference can be appreciated is 
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in cases where the systems at hand are subject to noise with some symmetries and we are 
interested in a fraction of the operators acting on the whole Hilbert space of the systems. For 
example, the polarization of photons moving in an optic fiber undergo collective rotation. 
When a symmetry of the noise exists, it is possible to find conserved quantities and thus 
invariant, or fixed, operators. These operators form the noise commutant A' of the noise. 
When the channel determined by the transmission is unital (i.e. the unit matrix is 

n 

preserved), the commutant is a C*-algebra IQ". Such an algebra has a block diagonal 
matrix form 

* 





where the compression to each block forms a full matrix algebra. But some of these blocks 
may be 'linked' in a sense that we make precise below. Once the block matrix form and the 
various links between blocks have been found, the explicit form of the algebra is revealed. 
In particular, there are unique positive integers rik, rrik > 1 for k = 1, . . . ,d such that the 
algebra is unitarily (spatially) equivalent to an orthogonal direct sum 

where is the full x Uk matrix algebra on n^-dimensional Hilbert space and Im^. 

is the identity operator on m^-dimensional space. A tensor product Ain^ ® in this 
decomposition corresponds to blocks of x matrices in (Q) which are linked, and 
is referred to as an ampliation of A^n^ when > 2. From the representation theory 
perspective, the integers correspond to the multiplicities of the identity representations 
of Airik which determine the structure of A'. 

If there are invariant operators which can be decomposed into projectors of rank one, 
so they are supported on blocks with = 1, then they correspond to invariant subspaces 
and the usual notion of qubits. Otherwise the invariant operators correspond to noiseless 
subsystems on higher dimensional spaces. This can be read off the structure of the com- 
mutant. For instance, in the above discussion there will be a copy of every operator from 
A4nk contained in the Aink ^ component of the decomposition. Thus the structure of 
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the commutant indicates how quantum information can be encoded to be preserved under 
the influence of the noise. 

In this paper we study how to find the structure of the commutant once the noise (Kraus) 
operators are known. A constructive proof of this structure is given in the case of unital 
evolution. The paper is structured as follows. In sections 2 and 3 we review background 
concepts. In particular, we recall basic properties of quantum channels and results from 
on the fixed point set of a unital channel. Section 4 contains the core theoretical component 
of the paper. Based on operator theory and operator algebra methods, we present a detailed 
algorithm for computing the noise commutant of a unital channel. We apply this method 
in Section 5 to find the commutant for a number of simple illustrative examples which are 



related to the phase damping channel 11 , 



Section 6 contains a short discussion on 



how to turn a non-unital channel into a unital one, when the noise operators satisfy certain 
constraints. In the final three sections we conduct an analysis of special cases of collective 



noise channels which arise from collective rotation 



Q Q Q Q, y , y, Q 



I . Specifically, 



section 7 includes an introductory discussion on this class, and in the final two sections we 



use the algorithm to explicitly compute the commutant for the 3-qubit [l^, 14,13, 
4-qubit channels from this class. 



22| and 



II. CHANNEL PRELIMINARIES 



Let 7i be a finite dimensional Hilbert space and let B{T-C) be the set of all bounded 
(continuous) operators acting on Ti. A quantum channel (or quantum evolution) on 7i is 
a linear map $ : B(H) — » BiH) which is completely positive and trace preserving. We 
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mention the texts |22l . I23l | for basic properties of completely positive maps, and 
introduction to channels in quantum information theory. 

It is well-known that to every completely positive map $ there corresponds a 

(non-unique) family of operators {Ai, . . . , An} in B{T-l) which determine the map $ through 
the equation 



$(T) = J] AkTAl for T G B{n). 



(2) 



fc=i 
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Trace preservation of $ is equivalent to these operators satisfying 

n 
k=l 

where 1 is the identity operator on Ti. (This is also equivalent to the 'dual map' of $ being 
unital.) We say $ is unital if also, 

n 

$(]i) = j]A,4 = ii. 

k=l 

For a quantum channel $, the family of operators {A{\ determining $ as in are called 
the Kraus noise operators for the channel. 
The fixed point set for $ will be denoted by 

Fix($) = {T G Bin) : $(r) = T}. 

We shall also let A be the so-called interaction algebra generated by the noise operators, 

A = k\g{Ai,...,An}. 

In other words, A is the set of all polynomials in the typically noncommuting variables 
Ai, . . . , An- Note that as an application of the Cayley-Hamilton theorem from linear algebra, 
the degree of any such polynomial may be reduced below some uniform bound. The noise 
commutant A! of A is the set of all operators which commute with every operator in A] it 
is the algebra 

a: = [T e B{n) : TA = AT for all A e A] 
= {TeB(n) : TAk = AkT for A;=l,...,n} 

Remark II. 1 The interaction algebra, and hence the noise commutant, are usually defined 
as f-algebras However, in the case considered here, that of unital noise, the discussion 
below shows that our 'nonself adjoint' definitions are actually a red herring, the algebras turn 
out to be selfadjoint. 

In this paper, we shall focus our attention on unital quantum channels. These are some- 
times called bistochastic channels in the literature. For such channels, it was proved in [igj ] 
that: 
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(i) ^ is a f-algebra; in other words, A belongs to A if and only if belongs to A. 
Furthermore, the algebra A is independent of the choice of noise operators which 
determine the channel as in (j2)); the same algebra is obtained whatever the choice of 
operators {Ak}. 

{ii) The commutant A! is a f-algebra which is equal to the fixed point set of $, 

A! = Fix($). 



These results are recalled in the next section. Taken together, they may be regarded 
as a structure theorem for the fixed point set of a unital quantum channel. In particular, 
the fixed point set Fix($) (which is just a f-closed subspace for general completely positive 
maps) is itself an algebra, hence closed under multiplication, equal to the noise commutant. 
But there is a tight decomposition theory for finite dimensional f-algebras (see section IIVI 
for details). This motivates the following fundamental problem which we consider in this 
paper: 

Problem II. 2 Given a unital quantum channel explicitly determine the algebra structure 
of the noise commutant Fix($) = A'. 

There are compelling reasons in quantum information theory for wishing to solve this 
problem. For instance, in the noiseless subsystem method of quantum error correction 

in : 

the structure of the noise commutant can be used to encode quantum 



a lid. 
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information which is immune to the noise of the corresponding channel. 



III. STRUCTURE THEOREM 

In jiol it was discovered that many properties of a channel are married to properties 
of the noise operators. Recall that an orthogonal projection P on Ti reduces an operator 
T in B{l-C) when both subspaces PTi and (1 — P)7i = P-^Ti are invariant for T; that is, 
T{P1-L) C PH and T{P^n) C P^H. This is equivalent to the identity TP = PT being 
satisfied, or T having a block diagonal matrix form T = [qq] with respect to the orthogonal 
spatial decomposition 7i = PTi. © P^Ti. Further recall the operator inequality P < Q for 
a pair of self-adjoint operators on H means that (-P^||0 ^ {Q^WO ^ E H. The key 
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point in the following result is that invariant subspaces for the noise operators of a unital 
channel are always reducing subspaces. 

Lemma III.l Let $ : i3(7i) BCH) be a unital completely positive map and let P be 

a projection on Ti. Consider the following conditions. 

{i) $(P) = p. 

{ii) The subspace P7i reduces Ai, . . . , An; 

AkP = PAk for l<k<n. 

(m) $(P) < P. 

(iv) The subspace PTi is invariant for Ai, . . . , A^; 

AkP = PAkP for l<k<n. 

(v) $(P) > P. 

(vi) The subspace PH. is invariant for y4|, . . . , A\^; 

PAk = PAkP for l<k<n. 

Then the pairs of conditions {{i), {ii)), {{Hi), {iv)), and {{v), {vi)) are each equivalent. More- 
over, if ^ is a unital channel then all six conditions are equivalent to each other. Further, 
given a projection P and a unital channel $, conditions {ii), {iv) and {vi) hold for one 
choice of noise operators for $ if and only if they hold for all choices. 

Let us outline the proof of {iv) =^ {ii) for a unital channel $, since this is the aspect of 
Lemma IIII.ll used in the proof of Theorem IIIL4I Suppose {iv) holds for a projection P in 
B{n). Then 

< $(P) = ^AiPAl = P$(P)P < P$(]1)P = PIP = P. 

i 

As < $(P) < P, trace preservation can be seen to imply through an operator theory 
argument that $(P) = P. More generally, if P is a positive contraction operator with 
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R < P for some projection P and R has trace equal to the trace of P, then in fact R = P. 
Thus $(P^) = $(1 — P) = P-*- since $ is unital, and for each j 

< {P^A]py{P^A]P) 

< {P^AlPy{P^AlP) = P$(P^)P = PP^P = 0. 

i 

Thus, P^^JP = for all j, and hence PH is a reducing subspace for ^41, . . . , 74„. 

We mention that this lemma gives a simple characterization of rank one projections fixed 
by a unital channel. 

Corollary III. 2 Let P = |^)(^| be a rank one projection in B{T-C). Then the following are 
equivalent: 

(i) $(P) = p. 

(ii) The range vector C, for P is a joint eigenvector for the noise operators Ai, . . . , An; 

^ = Afc ^ for some G C and k = 1, . . . ,n. 

Proof. Condition {ii) here is simply a restatement of {iv) from the lemma, in the special 
case of a rank one projection. ■ 

We recall some terminology from linear algebra before proceeding. 

Definition III. 3 Let T be a normal operator in i3(7i), in other words TT^ = T^T. Then a 
spectral projection for T is a projection onto an eigenspace for T; that is, it is a projection 
onto a subspace of the form E Ti. : T^ = X^} for some A G C. The Spectral Theorem from 
linear algebra states that every normal operator T decomposes as a sum T = A1P1+. . .H-A^Pr, 
where Aj G C are the distinct eigenvalues for T and Pj the corresponding spectral projections 
(which necessarily have mutually orthogonal ranges). 

Clearly the fixed point set Fix($) contains the noise commutant A'; if T belongs to A' 
then 

n n 

$(T) = A,t4 = t(5^ A,4) = Tm = T. 

k=l k=l 

The converse inclusion holds for unital channels. We mention that the following theorem 

n 

was partly motivated by work from [6] on infinite dimensional operator algebras and dilation 
theory. The proof in (l9| is somewhat technical, but after submission we discovered a simpler 
proof of this fact, and hence, for completeness, we shall present it here. 
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Theorem III. 4 jlylJ Let $ : B(T[) B{TC) be a unital quantum channel determined as in 
(0) by operators {Ai, . . . , An}. Then the fixed point set for ^ is a ^-algebra and is equal to 
the noise commutant, 

Fix{<^)=A'={A,,...,Any, 

Further, the interaction algebra A generated by Ai, . . . , An is a \-algebra which depends only 
on $ via the equation 

A = Fix($)'. 

Proof. Notice that the theorem will be proved if it is shown that A' = Fix($). Indeed, A' 
is an algebra by definition and Fix($) is f-closed since $ is a positive map. Moreover, if this 
identity holds, the algebra A depends only on $, not the choice of noise operators, since 
A = {A'Y = Fix($)'. (This is a special case of von Neumann's double commutant theorem 
from operator algebras.) 

As observed above, the commutant A' is contained in Fix($) because $ is unital. To see 
the converse, first observe that if the only fixed points are scalars, then Cll C C Fix($) = 
CI, and the two sets coincide. Otherwise, let T be a non-scalar operator in Fix($). Without 
loss of generality we may assume that T is self-adjoint; in other words, T = T"!". Indeed, 
$(rt) = <l>(T)t = Tt by positivity of $, and hence both the real. Re (T) = i(T + Tt), 
and imaginary, Im (T) = ^(7" — T"!"), parts of T = Re (T) + ilm (T) are fixed by $. Since 
T is non-scalar, at least one of these self-adjoint operators is non-scalar. Furthermore, by 
replacing T with T + | |T| |1, we may assume T > is a positive operator. Let < Ai < A2 < 
. . . < Ar be the eigenvalues for T, and Pi,...,Pr the corresponding spectral projections. 
Then T = AiPi + . . . + X^Pr by the Spectral Theorem, and < T < ||T||]1 = A,.]l. Let 
Hr = Pr'H = G : = A^^} be the eigenspace for the extremal eigenvalue A^. If ^ is a 
non-zero vector in Hr, observe that 

n 

\rm) = {Tm) = Y.^A,TAim) 

k=l 
k=l 

n 
k=l 

The only way this can happen is if each of the inequalities is actually an equality; in par- 
ticular, AI^ belongs to TCr for k = 1, ... ,n. Thus Hr = Pr'H is an invariant subspace for 
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A\ and so it is a reducing subspace for A by Lemma IlII.ll Therefore, it follows that Pr 
belongs to the commutant A', and hence to the fixed point set Fix($). Thus the self-adjoint 
operator T — X^Pr belongs to Fix($), and we may iterate this argument to find that each of 
Pi, . . . ,Pj. belongs to A'. Hence T belongs to A', as required. 

We have shown that every T = T"!" in Fix($) also belongs to A'. But Fix($) is a self- 
adjoint subspace, and hence spanned by its self-adjoint part. It follows that the fixed point 
set Fix($) is contained in the commutant A', and this completes the proof. ■ 



IV. COMPUTING THE COMMUTANT 

In this section, using operator theory and operator algebra techniques, we present a con- 
structive process for explicitly computing the noise commutant A' for a unital quantum 
channel. This will be accomplished by first identifying the f-algebra structure of the in- 
teraction algebra A. We begin with a short discussion on the decomposition theory for 
f- algebras. 

A fundamental result from the theory of operator algebras is that every f-algebra A 
is f-isomorphic (an isomorphism which preserves adjoints) to a unique orthogonal direct 
sum of full matrix algebras. Specifically, there are positive integers > 1 such that A is 
f-isomorphic to 

Ac:iMn,® ...®Mnd- (3) 

We mention the texts P, 0, |3] for basics of f-algebras, or equivalently, finite dimensional 
C*-algebras. Another important result is that every representation of such an algebra de- 
composes into a direct sum of multiples of the identity representation on each summand. As 
a consequence, up to unitary equivalence A is given by a unique direct sum of 'ampliated' 
full matrix algebras, 

^~ ®]i„j©...©(;w„, ®]i^j. (4) 

In other words, there is a uni tary operator U such that WAU is equal to the decomposition 
in dH). See also 



28L |29| for discussions of these decompositions in connection 



with quantum information processing. 
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While it is important to know the f-isomorphism class of ^ as in Q, we really need the 
structure of A spatially as in More to the point, we need to identify the precise spatial 
locations of the matrix blocks ® ^m^- The integers > 1 are the multiplicities of 
the identity representation appearing in each summand associated with the corresponding 
representation of the f-isomorphism class model in 0. The invariants nfc,mfc are part of 
what is known as the Wedderburn structure theory for such algebras. When some > 
2, we say Ain.k ® llm^ is an ampliation of A^n^. The projections in A which, under the 
unitary equivalence (jH), correspond to the identity operators on the blocks J^uk ® 
in the decomposition are called the minimal central projections for A. Each of these 
projections will be a sum of projections onto some, but perhaps not all, of the matrix blocks 
of the same size in the block diagonal form (Q). 

The projections onto the individual matrix blocks in will be minimal inside A', in 
other words minimal ^-reducing, but may not belong to the algebra A. This is precisely the 
subtlety we must deal with in this process; an identification of 'links' between these minimal 
projections. Thus, broadly speaking, we must first identify the maximal (here meaning 
that the projections sum to the identity operator) family of non-zero minimal ^-reducing 
projections with mutually orthogonal ranges (part 1 below); it is easy to see that such a 
family of projections is unique. Then we must detect any links between these projections 
(part II below) in the following sense. 

Definition IV. 1 Let {Pj} be the (unique) maximal family of non-zero mutually orthogonal 
minimal ^-reducing projections. A subset {Qk}kes of this family is linked in A if the 
following two conditions hold. 

(i) The projection Q = J2kes Qk belongs to A. 

{ii) If Si is a proper subset of iS, then the projection YlikeSi does not belong to A. 

Let {Pj : j G iSfejfc be the linked subsets of projections from {Pj}- It follows that each 
algebra A{^ "^j^s,, forms a subalgebra of A. Notice that, in terms of (ji}, the cardinality 
of the sets Sk are equal to the ampliation multiplicities m^. 

Once the f-algebra structure of A has been identified up to unitary equivalence, we can 
easily deduce the f-algebra structure of A' = Fix($). For instance, if $ is a channel and the 
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spatial structure of A is given by (jH), then 

Fix(<l>) =A'-Y1 ®(-^n, ® ImJ' - J2 ® 
fc k 

Furthermore, the algorithm presented below will allow us to explicitly identify the locations 
of the matrix blocks, possibly ampliated, which belong to the noise commutant. 

We shall work our way through the finer points of the process, then finish by summarizing 
the main points. For the rest of this section we let $ be a given unital quantum channel acting 
on B{T-C), with noise operators {Ai, . . . , An}, and interaction algebra A = Alg{y4i, . . . , An}. 
The following assumption gives us the starting point. 

Assumption: The noise commutant A' = Fix($) has been 
computed as a linear span. 

The algorithm will take us from knowing the vector space structure of A', to identifying the 
algebra structure of A, and hence of A'. There are two broad components to the process: 

Part I. Find the family of non-zero projections inside the noise commutant A' which 
have mutually orthogonal ranges, sum to the identity operator, and are each minimal 
inside A'. 

Part II. Given the family of projections {Pj} from Part I, find the minimal central 
projections for A by computing all links between the Pj. This will give the spatial 
f-algebra structure of A, and hence of Fix($) = A'. 

A. Part I 

As $ is unital. A' = Fix($) always contains the scalar operators. If A' = CI, then the 
singleton {11} is the only (non-zero) minimal central projection for A. Otherwise, let T = T'^ 
belong to A! \ CH. The spectral projections for T provide the first step toward the maximal 
family. 

Lemma IV. 2 The following are equivalent for an operator T = T'^ in B{T-C): 
(i) $(T) = T. 

{ii) Every spectral projection PofT satisfies $(-P) = P ■ 

12 



(m) Every spectral projection PofT belongs to A' . 

Proof. Condition (ii) and {iii) are equivalent by Theorem IIII.41 The equivalence of the 
first two conditions immediately follows from standard operator algebra for f-algebras. We 
provide a proof for the sake of brevity. The implication (ii) =^ (i) follows from the Spectral 
Theorem of linear algebra; T decomposes as a sum T = AiPi + . . . + A.,.Pr where Ai, . . . , A,, 
are the distinct eigenvalues for T, and Pi, . . . ,Pr the corresponding spectral projections. On 
the other hand, if $(T) = T, then T belongs to the commutant A' by Theorem 1111.41 But 
again by the Spectral Theorem, the spectral projections of T commute with everything that 
commutes with T. Hence (i) =^ (ii) follows since A' = Fix($). ■ 

Thus, if we let {Pj} be the (mutually orthogonal) spectral projections of T, then Pj 
belongs to Fix($) = A' for each j. Observe that for all j we may define a map $j : 
B{PjH) BiPjH) by 

n 

^j{S) = Y,A,,SAl, 

k=l 

where 

^k,j = AkPj = PjAk = PjAkPj. 

It should be kept in mind that this is a slight abuse of notation; we should really write Akj 
as the restriction Akj = Aklp^n- These are unital channels since 

n n 
k=l k=l 

A similar argument can be used to verify trace preservation for $j; that is, Ylt=i ^k = 
Ip^n- Proposition II V .41 below uses an analysis of these maps to give a test for deciding if Pj 
is minimal. 

Lemma IV. 3 Let \E' : i3(7i) —>■ B{l-L) he a unital quantum channel, with noise operators 
Ci, . . . , Cm, and let C = Alg{Ci, . . . , Cm}- The following are equivalent: 

(z) Fix(^) =C' = Cln- 

(ii) There are no non-trivial reducing subspaces for the algebra C. 
Proof. This immediately follows from Theorem 1111.41 since Fix(\l/) = C and C is a f-algebra. 
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Proposition IV. 4 Let \I/ : B{T-C) B{T-C) be a unital channel with noise operators {Ci\, 
and let P G i3(7i) he a projection which reduces the noise operator algebra C for Then 
the unital channel \l/p : B{PH) B{PH) given by \I/p(S') = J2i^i,P ^ ^1 P' where Ci^p = 
CiP = PCi, satisfies 

Fix(*p) = {CPy = PC P. 
Furthermore, P is a minimal projection in C if and only if 

Fix(^p) = PC'P = CP. 

Proof. Again, the statement (CP)' = PC'P is a notational convenience. We really mean 
that the compressed commutant, PC'|p7^, is equal to the commutant of the restriction, 
{C\pn)' = {PCIp-h)', inside B{Pn). To see this, first let T belong to C. Then 

{PTP){aP) = P{TQ)P = PiaT)P 

= {PCiP){PTP) = {aP){PTP), 

so that [CPy contains PC'P. Conversely, if T = PTP is in [CP)', then 

TQ = TPC, = T{QP) = {QP)T = QT, 

and T belongs to PC'P. Thus Theorem 1111.41 gives us 

Fix(^p) = {C,,p}' = {CiP}' = {CP)' = PCP 

The rest of the result follows from the previous lemma. ■ 

If we apply this result for P = Pj, we see that Pj is a minimal projection inside A' if 
and only if PjA'Pj = CPj. Thus, since we are assuming the noise commutant A' is known 
as a linear span, we can simply compute PjBPj for a set of vector space generators B for 
A' and check if each PjBPj belongs to CPj. If this is the case, then Pj is minimal and 
it belongs to the maximal family of projections we seek. On the other hand, if there is a 
Bq in A' such that PjB^Pj is not in CPj, then Pj is not minimal; however, in this case we 
may iterate the process outlined above. Indeed, $j is a unital channel on B{PjT-l) with non- 
scalar fixed points, in fact PjB^Pj G A'Pj = Fix($j) is such an operator. Thus, we may use 
PjBoPj to obtain a self-adjoint operator in Fix($j), and continue the process by examining 
the spectral projections of this operator, which are, of course, each a subprojection of Pj. 

14 



As Ti. is finite dimensional, this process will eventually terminate. We will be left with the 
maximal family (since they sum to the identity operator) of mutually orthogonal non-zero 
minimal projections in A', and this completes Part I. 



B. Part II 

We emphasize the importance of distinguishing the various links between projections in 
the maximal family from Part I by presenting the following simple example. 

Example IV. 5 The point is that we must distinguish between bonafide orthogonal direct 
sums on the one hand and ampliation algebras on the other by finding linked projections. 
Let 



j fa o\ 








\[o aj 





C = i \ I : a G C > ~ Clls ~ C ® lis, 
be a scalar ampliation algebra and let T) be the unlinked orthogonal direct sum given by 

V = {\ : a,6G C > ~C1©C1, 




where both of these f-algebras are diagonal with respect to the orthonormal basis {|0), |1)} 
for C^. 

The family {|0)(0|, |1)(1|} of rank one projections is the maximal family from Part I 
for both algebras. However, the commutant of the scalar ampliation algebra C = CUs is 
unitarily equivalent to C ~ M.2, whereas the commutant of the orthogonal direct sum 
T> ^ CI © CI is also V ~ CI © CI, the direct sum of two copies of the one dimensional 
scalar algebra. Notice that the projections |0)(0| and |1)(1| are linked in C, since C contains 
I2 = |0)(0| + |1)(1| and neither of these subprojections belongs to C. On the other hand, 
they are unlinked in V because both projections already belong to T>. 

Before continuing, let us recall a helpful identity from operator algebras. As a f-algebra, 
A is equal to its own double commutant, A = {A')' = A". This is a special case of von 
Neumann's double commutant theorem 0,0,|2^. Now, let {Pj} be the family of projections 
for A obtained in Part I. It follows from the discussion at the start of this section that links 
can only occur between Pj of the same rank. Hence for k > 1, let {Pj,k} be the projections 
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amongst {Pj} with dimPj a;?^ = k. Further let Sk be the index set for this set of projections. 
Suppose A' = span{i?i, . . . ,Br} with Bi in B{7i). 

bmce A = {AJ = (span{Bi,...,5,})' , we may check if projections of the form = 
Yljes ^hki where S C Sk-, belong to A simply by computing the commutators 

[Ps,Bi] = PsB,-B,Ps for z = l,...,r. 

The projections {Pj,k\jes are linked in A precisely when [Ps, Bi\ = for each i, and, there 
is no proper subset Si of S such that [P^^,!?^] = for each i.As Ti is finite dimensional, 
there are just finitely many computations required here. Thus, we will eventually exhaust 
all possibilities and discover the links between the projections {Pj,k}- 

Testing for links between the rank one projections in {Pj} is easier than the general 
case, hence we present it separately. Recall a rank one projection Q = belongs to 

A' = Fix($) if and only if ,^ is a joint eigenvector for Ai, . . . , An- 

Lemma IV. 6 Let {Qk} be the rank one projections from the set {Pj}. Let Qk = |^A:)(^fc| 
and suppose Xik are scalars with 

Ai^k = Kk^k for i = l,...,n. 

Define a function f : {Qk} — *■ C" by 

f{Qk) = (Alfc, . . . , Xnk)- 

If X in is such that /~"^({A}) is nonempty, then the projections in /~^({A}) are linked 
in A. 

Proof. Let /~^({A}) be a nonempty subset of {Qk}- Observe that if ^k and ^/ are range 
vectors for projections Qk and Qi in /^^({A}), then they are eigenvectors for A and we have 

{Aikm) = {Ami) for all A e A. (5) 

Indeed, f{Qk) = fiQi) and hence these inner products agree on a set of generators Ai, . . . , An 
for A. Suppose Q is a projection in A which is the sum of some, but not all, of the projections 
in /^^({A}). Then there is a pair of vectors ^k and ^i, each a range vector for a projection 
in /~^({A}), such that Q^k = (,k and QC,i = 0. In particular, the identity (0) does not hold 
for Q, and consequently Q does not belong to A. 
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For a nonempty set /~^({A}), let Qx be the sum of all projections inside /~^({A}). It 
remains to show that each Qx belongs to A. Prom the structure theory for f-algebras 
discussed at the start of this section, we know that for each m the projection Xlje^m 
belongs to where Sm is the index set for all projections in {Pj} of rank m. Let P = 
^^g^^ Pj = |Cfc)(Cfc| be the projection obtained for m = 1. Then it is easy to see that 
the compressions AiP — PAi, for i = 1, . . . , n, are mutually commuting operators inside 
A which are each diagonal with respect to the basis {^k} for the range of P. Hence it 
follows that AP is a subalgebra of A which, up to unitary equivalence, has a block matrix 
decomposition of the form 



AP 







y an^KdJ 

Further, under this unitary equivalence, each projection onto a block l^j corresponds to 
a projection which is the sum of projections from a particular nonempty /~^({A}). Since 
complete freedom is allowed for the scalars a„ ., we see that every Qx belongs to A. ■ 



Remark IV. 7 As a prelude to some of the examples below, we note that the case of an 
abelian algebra A corresponds to the case that the family of minimal ^-reducing projections 
are all rank one. 



C. Summary of the Algorithm 

We present a number of applications in the rest of the paper, but let us finish this section 
with a brief technical summary of the algorithm. Recall we have assumed that a set of vector 
space generators is known for Fix($) = A' = span{i?i, . . . , Bj.}. 

Part I. 

(i) If A' ^ CI, then choose T = Tt inside A' \ CI. 
(ii) Compute the spectral projections Pi, . . . ,Ps for T. 

(Hi) For each P ^ Pj, compute PB^P, . . .,PBrP. If every PBiP belongs to CP, then 
PA'P — CP and P is a (non-zero) minimal reducing projection for A, and it belongs 
to the maximal family. 
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(iv) If some B = Bi is such that PBP does not belong to CP, then P is not minimal 
^-reducing and the compression channel has non-scalar fixed points; PBP for 
instance. Then proceed through steps (i), (ii) and (iii), perhaps a number of times, 
to find minimal ^-reducing mutually orthogonal projections supported on P. 

Part II. 

(i) Given the family {Pj} from Part I, group together those projections with the same 
rank. 

{ii) For the projections {Pj,k}jeSk '^^ rank k > 2, test for links between the Pj^k by evaluat- 
ing the commutators [Ps, Bi] where Ps = J2j(^s with S C Sk- Then {Pj^k : k E S} 
are linked in A if and only if [P^, Bi\ = for each i, and there is no proper subset of S 
with this property. The projections of rank k = 1 may be tested using Lemma fiV. 61 
The projections Ps defined by linked sets of projections will form the family of minimal 
central projections for A. 

(Hi) The various links will give the f-algebra structure of A up to unitary equivalence, and 
hence of A' = Fix($), along with the precise locations of matrix blocks inside A'. 



V. PHASE DAMPING AND RELATED EXAMPLES 



In specific cases the computations involved with the method described above can be 
less cumbersome than what is suggested. We begin with a number of simple illustrative 
applications of the process. In each of these cases the commutant is already well-known, or 



at 



11 



east easy to find, but we shall demonstrate how to compute it using the algorithm. See 



IJ] for related examples. 



Recall the Pauli matrices are given by 



Is 







We regard these as matrix representations of operators acting on with respect to an 
orthonormal basis {|0), |1)} corresponding to the base states of a given two- level system. 
In the case of higher dimensions, for instance with basis {|00), |01), |10), |11)}, we use 
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notation such as Z\ and Z2 to denote, respectively, the tensor products oi Z®\2 and \2®Z 
acting on C^. 



Example V.l Let p be a positive real number, < p < 1. Let Ai, A2 be operators on 
defined on the standard basis by: 

= (v^r^)]i2 = 



and 

A2 = {^)Z = ^ 




1 
-1 



Clearly ^Yl!k=\ ^fc^I = ^2 = Y^k=\ ^I^fc? ^o {Ai, ^2} are noise operators for a unital channel 
$ acting on i3(C^). The quantum operation correspondin g to this channel is equivalent to 
the phase flip or phase damping operation on single qubits j2l|. Its efi^ect is to flip the phase 
of |0) relative to that of and vice-versa, with probability p. 
A simple computation shows that 



Fix(<l>)=^' = {^1,^2}' 
a 



a,6G C > ~C1© CI 



b 

= A\g{Z} = A. 

Let us see how this structure arises through the algorithm. First choose a non-scalar operator 
in A'] in this case Z will suffice. The rank one projections Pq = |0)(0| and Pi = |1)(1| are 
the spectral projections for Z, hence they are minimal inside A' since the vectors |0) and 
|1) are joint eigenvectors for Ai,A2. Thus the family {Po,Pi} forms the maximal family 
(12 = Pq + Pi) of non-zero mutually orthogonal minimal projections inside the commutant 
A'. 

Further observe that 

Ai\0) = VT^\0) j Ai\l) = \l) 

and < 

A,\0) = ^\0) [A2\1) = -^\1) 

The 2-tuples {^/T^^, ^/p) and (-^/l — p, —y/p) are distinct, so the projections Pq and Pi 
are not linked by Lemma flV.GI Hence we have A = APq © APi, with APq and APi both 
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subalgebras of A each unitarily equivalent to the one dimensional scalar algebra Cl. Thus, 
the commutant A' — (APq)' + (APi)' — PqA'Po + PiA'Pi is also unitarily equivalent to 
Cl © Cl, and the form of Fix($) = A' given above is evident. 

Example V.2 Let p be a positive real number, < p < 1. Let Ai,A2 be operators on C^ 
defined on the standard orthonormal basis by: 

/l o\ 
10 



Ai = y/l-p {h h) = V^-P 



10 
\0 1/ 



and 



/l o\ 
0-100 
0-10 
\0 1/ 

Once again, it is clear that these noise operators determine a unital channel $ on B{C^) 
We may compute 

Fix($)=^' = {A^,A2}' 

= A\g{Z^,Z2,X^X2 + Y,Y2} 
/oii ai4\ 

022 023 

032 O33 

^ \a4i 044/ 

By observation, we can see that A' is unitarily equivalent to the orthogonal direct sum 
A' ~ M2 © A^2- To see this through the algorithm, first choose a non-scalar self-adjoint 
operator in Fix($) = A'; 

/l \ 

10 
0-10 
\0 -1/ 



: aij G C 



> . 



Zi = Z (g) I2 
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is such an operator. The spectral projections for Zi, corresponding to eigenvalues A = 1 and 
A = — 1, are given by 

Pi = |00)(00| + |01)(01| and P_i = |10)(10| + |11)(11|. 

Let Pkj — \kj){kj\ for k,j E {0, 1}. The projections Pi and P_i belong to A', but in this 
case they are not minimal inside A'. Indeed, observe that 

P.A'Pi = {APi)' ^{AiPi,A2Piy 
' /an 0\ 

022 


\ 0/ 

and similarly for P_i. Following the algorithm, we let 

/l o\ 
0-100 


yo 0/ 

be a non-scalar operator in Pi^'Pi. The spectral projections for this operator are 



: an, 022 G C 



Z1Z2P1 — P1Z1Z2 



00 



|00)(00| and Pqi = |01)(01|. 



These rank one projections are minimal inside A'. Similarly, we find the subprojections 
Pio = 1 10) (10 1 and Pn = |11)(11| of P-i are minimal inside A', and hence the maximal 
family of non-zero mutually orthogonal minimal projections in A' is given by 

It remains to check for links between these minimal projections. Observe that 



Ai|00) = |00) 

A2I00) = VP |oo) 



A|01) = VT^IOI) 

A2|oi) = -Vp|oi) 



Ai|io) = vr^lio) 

A2|10) = -Vp|10) 



= 111) 
A2|ii) = Vp|ii) 
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Thus, upon comparing the corresponding 2-tuples — p, ^/p) and (^^1 — p, — a/p), we see 
that the pairs {Pqq, Pu} and {Pqij-Pio} are hnked in A by Lemma flV.61 In particular, 
Poo + Pn and Pqi + ^'lo belong to ^ = A{Poo + Pii) © A{Poi + Pio) and A is unitarily 
equivalent to ^ ~ C112 © C112. 

Hence, Fix($) = A' ^ A^2 ©A^2, and the locations of the M2 blocks inside Fix($) = A' 
are explicitly given by the algorithm. Namely, one block of 2 x 2 matrices inside A' = Fix(<l>) 
is given by all T in i3(C^) such that T = (Pqo + Pii)T(Poo + Pn), and similarly the other 
M.2 block is supported on Pqi + Piq. The blocks are disjoint in A' in the sense that if T 
belongs to A^2, then an operator T = (Pqo + Pii)T(Poo + Pn) can be defined in A', and 
similarly for the Pqi + Pio block. 



Example V.3 Let p be a positive real number, < p < 1. Let Ai, A2, A3, A4 be operators 
on defined on the standard basis by: 

/l o\ 

10 



A, = {l-p) (]l2 © 112) = (1 - p) 



A2 = (Vp(i-p)) z, = v/p(i-p) 



A3 = (Vp(i-p)) Z2 = Vp(i-p) 



= (p) Z1Z2 = {p)Z ®Z = p 



10 

\o ly 

/l \ 
10 
0-10 
\0 -I J 

/l \ 

0-100 
10 

\o -ly 

/l o\ 

0-100 
0-10 

Vo ly 
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These noise operators determine a unital channel $ on B(C^). The fixed point algebra may 
be computed as 

Fix($) = 



= < 



{^1,^2,^3,^4}' 

Va 0\ 
6 
c 



: a,b,c,d E 



dj 

~ cieciecieci 

Running through the algorithm, we first choose Zi — Z <S) ^2 — Zl as a non-scalar 
operator in A'. As in the previous example, the spectral projections for Zi are given by Pi 
and P-i. Once again, neither of these projections is minimal inside A'. For instance, for 

= |00)(00| + |01)(01| the fixed point set of the restricted channel ^p^{S) = Y.t=i ^A^A,!^ 
where 74^^! = AkPi, satisfies 

Fix($pJ = {AP,y = P.A'Pi 
' /a 0\ 
6 


\o oy 

and hence contains non-scalar operators. 

It follows that the projections {Pqo, Pqi, Pio, Pu} again form the maximal family of non- 
zero mutually orthogonal minimal projections in Fix($) = A'. But unlike the previous 
example, there are no links between these rank one projections. Indeed, observe that the 
numeric 4-tuples given by the equations 



= < 



a,6 e C 



' Ai\00) = 


--{l-p) ICQ) 


Ai\01) -- 


= (l-p) |01) 


^|00) = 


= Vp(i-p)|oo) 


A^lOl)-- 


= Vp(i-p) |oi) 


AalOO) = 


Vp(i-p) |oo) 


AslOl) -- 


= |oi) 


^ A4IOO) = 


.p|00) 


A^\01) -- 


= -P |oi) 


' Ai\10) = 


|10) 


' Allll) 


= {l-p) 111) 


A2|10) = 


-Vp(i-p) |io) 


A2\n) 


^-^/p{l-p) 111 


A3IIO) = 


v/p(i-p) |io) 


Asln) 


= -a/p(i-p) 111 


^ ^4|10) = 


-P |io) 


. ^4|11) 


=plll) 
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are pairwise distinct. Thus each projection Pij belongs to ^ = APqq + ^Poi + APio + APu, 
and A is unitarily equivalent to the unlinked orthogonal direct sum A ~ C1©C1©C1©C1. 
Hence, Fix(<l') = ^' ~ CI © CI © CI © CI as well. 

VI. UNITIZATION OF CHANNELS 

Before proceeding we present a discussion which will be useful in the sequel. Recall 
from the Spectral Theorem that every normal operator T in B{T-l) decomposes as a sum 
T = AiPi + . . . + XrPr where Aj are the distinct eigenvalues for T and Pi the corresponding 
spectral projections. The functional calculus for T is determined by all complex-valued 
functions / which are defined on the eigenvalues of T. For each such function, there is an 
operator /(T) in the functional calculus for T defined by /(T) = /(Ai)Pi + . . . + /(Ar)-Pr- 

Lemma VI. 1 Let T be a normal operator, and suppose f is a function such that f{T) 
belongs to the functional calculus for T . If f is injective on the set of eigenvalues for T , 
then the commutant {T}' coincides with the commutant {/(T)}'. 

Proof. Suppose T = AiPi + . . . + XrPr is the spectral decomposition of T. Then, by 
the Spectral Theorem, the commutant of T may also be realized as the commutant of 
its spectral projections; {T}' = {Pi, . . . ,Pr}'. By hypothesis, /(Aj) = f{Xj) implies i = j. 
Thus, it follows that the spectral projections for f(T) = /(Ai)Pi + . . . + f{\r)Pr are precisely 
{Pi, ... , Pr}. Hence, {T}' = {/(T)}', as claimed. ■ 

As an immediate consequence of Theorem lIII.4t we obtain the following result for general 
quantum channels. 

Corollary VI.2 // {Ai, . . . ,An} are normal operators, and if, for each k, f{Ak) is in the 
functional calculus of Ak and is injective on the set of eigenvalues for Ak, then 

{Ai,...,A4' = {/(Ai),...,/(A„)}'. 

In particular, if 

n 

<l>(X) = 5^/(Afc)X/(Afe)t 

k=l 

defines a unital quantum channel, then 

Fix($) = {Ai,...,A„}'. 
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Hence this gives a way of 'unitizing' a quantum channel, when the noise operators satisfy 
the hypotheses of Corollary IVI.21 in such a way that the noise commutant of the new unital 
channel is completely determined by the commutant of the original noise operators. We 
shall make use of this lemma in the following section with the function / given by a multiple 
of the natural exponential function. 



VII. COLLECTIVE NOISE CHANNELS 



We next consider a class of examples which are more involved technical 

iterature 



I: 



The^qubit 
and 



example from this class has been analyzed previously in the literature jiol . 14 , l^l 
also experimentally realized via nuclear magnetic resonance j27|. In particular, the algorithm 
allows us to explicitly identify the structure of the noise commutant for this channel. We 
will also conduct this analysis on the 4-qubit case, but first we outline the general n-qubit 
case. 

Let n > 3 be a fixed positive integer. Let X be defined on C^" by 

X = Xi + X2 + ...+X„, 

where 

Xi = X ® (nf "-^)) , X2 = ]i2 ® X ® (]if , . . . 

Similarly define Y = Yi + . . . + Yn and Z = Zi + . . . + Zn- Observe that all of these 
operators are self-adjoint. Let exp(-) be the complex-valued natural exponential function. 
This function belongs to the functional calculus of each of X, Y, Z, and hence we may define 
operators Ex,Ey, on C^" by 

E^ = ^ exp(zX), Ey = ^ exp(zY), E, = ^ exp(zZ). 

A unital channel $, called the n-qubit collective noise (rotation) channel, is then defined on 
BiCn by 

$(T) = E.TEl + EyTEl + E^TEl (6) 

As before, let A = A[g{Ex, Ey, Ez] be the algebra generated by the noise operators. The 
following result provides a useful computational device for this class of channels. 
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Proposition VII. 1 //$ is given by (0), then 



Fix(<l>)=^' = {E,,E„EJ' = {X,Y,Z}' 

= {X,Y}' = {X,Z}' = {Y,Z}'. 



Proof. The identities in the first line of this equation are a consequence of Theorem 1111.41 
and Corollary IVI.2I for the equality {E^^, Ey, E^}' = {X, Y, Z}' we use the fact that exp(i-) 
is injective on the set of eigenvalues for each of X, Y, Z. This can be seen by noting that 
the spectrum of Z consists of integers and that X, Y, Z are unitarily equivalent. 

The equation {X, Y, Z}' = {X, Y}' follows from the anti-commutation relations for the 
Pauli matrices since 



Remark VII. 2 Notice that we could also define a channel $ for n = 1 and n = 2. But 
simple computations in both cases show that Fix($) = {X, Y, Z}' contains no bonafide 
matrix blocks; Fix($) = C112 for n = 1, and Fix($) ~ Cl^ © CI for n = 2. Hence these two 
channels are not of interest from a quantum computing perspective. On the other hand, 
there is rich structure in the commutant for n > 3. 

VIII. THE 3-QUBIT CASE 

In this section we consider the n = 3 case for the collective noise channels. The matrix 
representations we use in this section are with respect to the standard ordered orthonormal 
basis for C*^ associated with the usual tensor product notation A^B = {aijB)ij for matrices. 
This is the basis 



We show explicitly that the algebra A is unitarily equivalent to the unlinked orthogonal 
direct sum of the full matrix algebra Ai^ together with the ampliation algebra A^2 ® ^2', 




The other equalities are similar. 



(1000), |001), |010), |011), |100), |101), |110), |111)}. 
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Thus the noise commutant contains an amphated copy of Ai2', 



Fix($) = ^' ~ {M2 h)' © (-M4)' 
~ (I2 ® M2) © C114 



A 







b 



■.AeM2,beC 



> 



b 







6 



v 



v 



We will verify this by using the algorithm to produce the maximal family of non-zero minimal 
^-reducing projections {P,Q,R}, and show that P,Q are linked inside A with rankP = 
2 = rankQ and ranki? = 4. In particular, this will imply that A = A{P + Q) (B AR with 

A^2 © ]l2 ^ A{P + Q) ^A and Mi^AR^ A, 

and thus Fix($) = A' = {A{P + Q))' ® {AR)' with 

M2 (S) I2 - {A{P + Q))' C A' and CI4 - {AR)' C A' . 

We will also explicitly identify a copy of the Pauli matrices sitting inside the ampliation 
A<2 © contained in A' = Fix($). 

We begin by computing A' as a linear span. By Proposition IVII. 11 we have A' = {X, Z}'. 
If we let B = {bij)i<i,j<8 and consider the commutators [B, Z] and [B, X], we find a system 
of linear equations determined by the bij which is satisfied precisely when B belongs to A'. 
Specifically, the system includes 30 equations with a total of 20 unknowns. Thus we may 
apply basic linear algebra (we use Matlab to put the corresponding matrix in its reduced 
row echelon form) to obtain a set of vector space generators for A' . As it turns out, there 
are five parameters which yield matrices Br, Bg, Bt, Bu, B^. That is. 



But B^ = B^. Hence it is clear from Parts I and II that, for the purposes of the algorithm, 
we may restrict our attention to the set {Br, Bs, Bt, Bu}- The matrix representations of 



A' = span{Br, Bs, Bt, Bu, By}. 
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these operators with respect to the standard basis for are given below. 



Bs = 



Bt = B] 
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1 1 
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Bu — 



/OOO 00 o\ 

1 0-1000 

000 000 

1 0-100 

0-10 1 000 

000 000 

-1 1 

\ooo oooy 

For Part I of the algorithm we shall conduct the spectral analysis on i?^ = Bl. The 
characteristic polynomial of Br may be computed as det(A]l8 — Br) — A^(A — 1)^(A + 2)^. 
Let Po,Pi,P-2 be the projections onto the eigenspaces for the eigenvalues A = 0, A = 1, 
and A = — 2. These projections belong to A' with rankPo — 2 — rankP_2 and rank Pi = 4. 
Define vectors {^o, ^o, C-2, ^7-2} in as follows: 



/o\ 

1 

-1 








1 



/o\ 







1 

-1 



C-2 = ^ 



1 

\7t 



/o\ 

1 
1 



-2 





/o\ 



2 


-1 
-1 

VO/ 



We note that these vectors are essentially the 'singlet-triplet' basis for the noiseless repre- 
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sentation as given in 



26j | and also used in [?, 10, Q. Then Pq and P_2 are given by 

Po = \^o){U + \Vo){Vo\ 



and 



P_ 





/O 








o^ 






1/2 -1/2 















-1/2 1/2 


















































1/2 


-1/2 









-1/2 


1/2 









^0 








V 




^-2){^^2\ + \V-2){V-2\ 










/o 














o\ 





1/6 1/6 - 


-1/3 














1/6 1/6 - 


-1/3 














2/3 


- 


-1/3 


-1/3 





-1/3 -1/3 


2/3 














-1/3 





1/6 


1/6 





-1/3 





1/6 


1/6 
















V 



Furthermore, Pi is computed as 

Pi 



Hg - Po - P_2 

/lo oo\ 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 

1/3 1/3 1/3 

\o oiy 
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Next we must examine the compressions of the commutant by these projections. First 
observe that P\Br = BjPx = XPx for A = —2, 0, 1 as the Pa are spectral projections for Br. 
We further compute 

P-2BsP^2 = P-2BtP-2 = P-2BuP-2 = P-2, 
PlB.P^ = P^BtPl = PAPl = 0, 

PoB.Po = PoBtPo = Po and PoB^Po = (^-0 Pq- 

Hence it follows that PxA'P\ = CPa for A = —2,0,1. Therefore, the minimal ^-reducing 
projections {P„2,Po)Pi} form the desired set of projections. 

It remains to check for links. Since linked projections have the same rank, the only 
possible link between these projections is P_2 + Pq. But notice that PqPu 7^ B^Pq, and 
hence 

Po ^ {A')' = A" = A. 

Thus, it follows that Pq and P_2 must be linked inside A. (One may check that Pq + P_2 G 
A = A".) This shows that A, and hence A', has the form claimed above. 

Finally, we shall exhibit a copy of the Pauli matrices inside A^2 ® I2 — (Po + P-2)'^'(Po + 
P-2) ^ A' = Fix($). We may accomplish this by first identifying a set of 2 x 2 matrix units 
{Ell, E12, E21, E22} inside (Pq + P_2)-^'(Po + P-2)- For this the algebra structure of A and 
A' may be used together with a spectral analysis of the operators {X, Y, Z}. Specifically, 
with the orthonormal basis {^0,^0,^-2,^-2} for PqC® © P-2C^ we may define 

Eii = \^o){^ + \Vo){^o\=Po 

Ei2 = \^0){^-2\ + \m){V-2\ 

E21 = 1^-2) (Col + 1^-2) (^?o| 

E22 = |C-2)(C-2| + \V-2){V-2\ = P-2 

Then the matrix representations of operators inside f Po + P_2)'^'(Po + P-2) with respect to 



the basis {^o, C-2, "/^o, ''7-2} (a 'canonical shuffle' (22, |2^ of the above basis) have the form 



(Po + P-2)^'(Po + P-2) ^ <( I ^ ° I : A e A^2 
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and the Eij are the natural matrix units in this decomposition. Thus, we may obtain versions 
of the Pauh matrices inside A' — Fix($) by defining 

/oooooooo\ 

1 0-10 
-1 1 
1-10 

\/30-ll 00 
1 0-100 
-1 1 

\ooooooooy 



X — Elo + Eol — 



Y = -iEio + iEoi = 



Ell — 
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One may verify directly that the anti-commutation relations are satisfied by {X, Y, Z}, and 
that this triple belongs to Fix($) ^ A' ^ {X, Y, Z}'. 
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IX. THE 4-QUBIT CASE 



In this section we determine the structure of A' = Fix($) for the collective noise channel 
with n = 4. We will show that A is unitarily equivalent to the unlinked orthogonal direct 



sum 



A-Ch® {Ms (g) h) 



M5, 



and hence 



Fix($) =A'-M2® {M3 ® h) © CI5. 



We can recognize a noiseless (decoherence-free) subspace |8|, |20|, |30| which supports the first 
summand in the decomposition of A and a noiseless subsystem in the second summand. We 
establish the above forms of A and A' by identifying the maximal family of projections from 
Part I of the algorithm, {P-3, P-i, Po-1/3, Po,i, Pi, P3}, and show that: P_3 Pi are linked in 
A with rankPfc = 1 for = —3, 1; the triple P_i, Po,_i/3, P3 is linked in A with the rank of 
each equal to three; and rankPo^i = 5. 

As for the 3-qubit case we may use elementary linear algebra to identify a set of vector 
space generators for A'. It turns out that it is sufficient in this 4-qubit case to obtain a set 
of self-adjoint generators which have real matrix entries. We find that A' is spanned by 10 
such operators. For succinctness we shall only display the generators used in the spectral 
analysis. Let 



Bn = Bn 



0000000000000000\ 

ooooooooooooooooN 

00001000 -1 0000000 
0000001000 -1 00000 
00100000 -1 0000000 
000000100000 -1 000 
000101000 -2 000000 
000 2 -1 -1 00 
00 -1 -1 2 000 
000000 -2 000101000 
000 -1 000001000000 
0000000 -1 00000100 
00000 -1 0001000000 
0000000 -1 00010000 
0000000000000000/ 

0/ 



\ooooooooooooooo 
This operator belongs to the commutant with characteristic polynomial det(A]li6 — Pq) = 

(A + 3)(A + l)3AS(A-l)(A-3)^ Consider the set of spectral projections {P_3, P_i, Pq, Pi, P3} 

corresponding to the eigenvalues Ai = —3, A2 = —1, A3 = 0, A4 = 1, A5 = 3. These 

projections belong to A' with rankP_3 = 1 = rank Pi, rankP_i = 3 = rankP3, and 

rankPo = 8. Let 



1 

2V3 



000101-200-2101000]* 
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,^1 = -[00010 -1 0000 -1 0100 o]* 
^-l^l = -^[OOOOOOOOOOO -10100 



e-1,2 



V2 

1 



000 -1 010000 -1 01000 



3 — -ij[00-l 0100000000000 ]* 

= --i=[oo-i 0-1 00020000000]* 
{3,2 = [000-10-1-2002101000]* 

^3 3 — -ig[0000000 -2 00010100]*. 

Then these projections may be computed as 



^-3 = iei)(6i 
^^1 = 1^4) (61 

P-i = I6,i)(6,il + 16,2) (6,2! + 16,3) (6,3l 

P3 = |6,l)(e5,l| + \^5,2){^5,2\ + |e5,3)(e5,3| 
Pq — I16 ~ P-3 — P-1 — Pi — P3 



As rank one projections, P-3 and Pi are minimal in A'. Further, compressing the genera- 
tors of A' to P_i and P3 reveals that PA'P = CP for both P = P_i and P = Pg- Hence P_i 
and P3 belong to the maximal family. On the other hand, the projection Pq is not minimal 
^-reducing. There are a number of ways to see this, including verifying that Pq^'Po 7^ CPq 
directly, but perhaps the easiest way is to observe that there arc vector space generators 
for A' which have all their eigenvalue multiplicities strictly less than 8 = rankPo. We could 
find the (unique) maximal family of projections by conducting a spectral analysis on any 
one of the non-scalar vector space generators of A'. Only two of the ten generators B for 
A' satisfy PqPPo 7^ 0. They are: 



Pi = Bl 
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1 






-1 
100 
0-20 

200 




1 





1 
000 




000 
000 

1 
000 

1 
200 
001 
2000 
1 -1 
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000 
000 
000 
000 
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-1000 
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000 
000 
000 
000 
001 
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and 



Bo 





1 









1 1 



0-10 

1 1 

-1 -1 



-1 -1 





.0 




0-1 



10 



1 
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-2 
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1 









10 

0-10 










-1 -1 



-1 -1 

10 1 





10 1 















0\ 
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0-10 













1 





° / 

0/ 



Following the algorithm, we shall consider the spectral projections for PoBiPq. We 
compute det(A]li6 - PqBiPq) = (A + l/SfX^X - if. Let Po,-i/3, ^o,o, A,i be the pro- 
jections onto the eigenspaces for Ai = —1/3, A2 = 0, A3 = 1. Since P^j- < Pqq and 
rankP(^ = 8 = rankPo,0) we have Po,o — Po~j we need only focus on Po,-i/3 and Po,i- 
Observe that rankPo,_i/3 = 3 and rankPo^i = 5. Let 

1 r It 
TJi — -= [00000001000101 -3 Oj 

2v3 



r ^t 

772 = — p [000 -10-1 100 -iioioooj 
1 r It 

773 — — — [o-aioioooioooooooj 
2v3 

Then Po,-i/3 and Po,i may be computed as 

^0,-1/3 = |^l)(^l| + |^2)(^?2| + \V3){V3\ 
Po,l = Pq — -fo,-l/3- 

For the generators B of A' other than Bi , B2 we have 

Po,_i/3(P)Po,-l/3 = Po,-l/3{PoBPo)Po,-l/3 = 0, 

and similarly Po,iPPo,i = 0. Moreover, 

1 4 

-fo,-l/3(-Bl)Po,-l/3 = — g-fo,-l/3 -fo,-l/3(-B2)Po,-l/3 = g-fo,-l/3 



and 



Po,l{Bl)Po,l = ^0,1 ^^0,1(52)^0,1 = 0. 



Thus PA'P = CP for both P = Po,-i/3 and P = Po,i, and hence Po,-i/3 and Po,i are 
minimal ^-reducing projections. 
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As lie = -P-3 + P-i + -Po,-i/3 + Po,i + -Pi + P3, the set of minimal ^-reducing projections 
{P_3, P-i, Po,-i/3, Po,i, Pi, P3} does indeed form the maximal family of projections we seek. 
It remains to find links. The commutators [P-3, B] for generators B of A' are not all zero. 

Hence, as the only rank one projections in this family, we may deduce that P_3 and Pi are 
linked inside A. Thus, P_3 + Pi G {AJ = A and the block ^(P_3+Pi) = (P_3+Pi)^ ~ CI2 
so that 

^'(P_3 + Pi) = (P_3 + Pl)^' ~ M2. 

Further, a check of commutators [P,B] for generators B of A' with P = P_i, P3, Po_i/3 
shows that none of these projections belongs to {A')' = A. It follows that P-i, P3, Po,-i/3 
are linked inside A, and thus 

A{P-i + P3 + Po,-i/3) ^ Ms^h 

~ 1^0 M3 

~ (Ms^h)' -A'{P-l + P3 + Po,-l/3)■ 
Fmally, we have APo,i ^ M5 and A'Po,i ^ (AI5)' = and the structure of Fix($) = A' 
is now apparent. 

We finish by identifying copies of the Pauli matrices belonging to the matrix blocks inside 
A' = Fix($). For the block 7W2 — ^'(P-3 + Pi), we may define matrix units in A' by 

En = P-3=\U){^~3\ Eu = \C-3){Ci\ 

E21 = \Cl){i-3\ E22=Pl=\il){il\ 

In particular, we may define the following operators inside ^'(P_3 + Pi) Q A' — Fix($): 

X — E12 + E21, Y — —iEi2 + i-E'21, Z — Ell — E22- 

We may also obtain matrix units inside A' for the copy oi Ai^^i^ — -4'(P_i + P3 + Po,_i/3) C 
A'. For instance, a set of 2 x 2 matrix units inside this block is given by: 

Eu = P-l = Yll=l IC-l,fc)(C-l,A:| E12 = Ylt=l \^-l,k) {^3,k\ 
E21 = Yfk=l \i3,k) {C-l,k\ F22 = P3 = Z)fe=l IC3,A;)(C3,A;| 

Thus, a second copy of the Pauli matrices inside A' — Fix($) is realized by defining 
X = P12 + P21 , -iFi2 + iF2i , Z^Fii- F22. 
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X. CONCLUSION 



Given a unital quantum channel $, we have derived a constructive proof for computing 
the exphcit algebra structure, as in (@)), of the associated interaction algebra A and noise 
commutant A'. We mention that Zarikian jsij has recently written Matlab algorithms 
which also utilizes much of the theory discussed here to compute, among other things, the 
structures of such operator algebras. An important subtlety in this process involves the 
detection of ampliations within these algebras which arise from linked minimal projections. 
In the case of non-unital noise the algorithm may be easily adapted to find the structures of 
A and A', since both are f-algebras by definition, but the connection with the fixed point 
set Fix($) = A' is no longer valid. 

As illustrations we worked through the process for several simple channels and the n = 3 
and n = 4 qubit cases of the collective noise channels which arise from collective rotations. 
We have also used our method to compute the noise commutant for the 5-qubit case, it is 
given by A' ~ Ole © (Ma ^ I4) ® {M5 ^ h)- 

While the algorithm is capable of analyzing any situation, in particular cases a more 
specialized and conceptual investigation may be necessary. For instance, when n >> the 
computations in the algorithm required to compute the noise commutant of the n-qubit 
collective noise channel make the problem infeasible. In forthcoming work, the authors and 
Pouhn [3| use such a conceptual investigation to derive the explicit structure of the noise 
commutant for the general n-qubit case of the collective noise channels arising from collective 
rotation. 
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